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ABSTRACT 



o 

■ Context. HD 181327 is a young Main Sequence F5/F6 V star belonging to the ft Pictoris moving group (age ~ 12 Myr). It harbors 
1^-^ ' an optically thin belt of circumstellar material at radius ~90 AU, presumed to result from collisions in a population of unseen 
planetesimals. 

Aims. We aim to study the dust properties in the belt in great details, and to constrain the gas-to-dust ratio. 
Methods. We obtained far-IR photometric observations of HD 181327 with the PACS instrument onboard the Herschel Space 
0^ ' Observatory*, complemented by new 3.2 mm observations carried with the ATCA** array. The geometry of the belt is constrained 
| with newly reduced HST/NICMOS scattered light images that allow the degeneracy between the disk geometry and the dust 
• , properties to be broken. We then use the radiative transfer code GRaTer to compute a large grid of models, and we identify the 
i-C . grain models that best reproduce the Spectral Energy Distribution (SED) through a Bayesian analysis. We attempt to detect the 
Oh- oxygen and ionized carbon fine-structure lines with Herschel/PACS spectroscopy, providing observables to our photochemical code 
ProDiMo. 

Results. The HST observations confirm that the dust is confined in a narrow belt. The continuum is detected with Herschel/PACS 
completing nicely the SED in the far-infrared. The disk is marginally resolved with both PACS and ATCA. A medium integration 
of the gas spectral lines only provides upper limits on the [OI] and [CII] line fluxes. We show that the HD 181327 dust disk consists 
of micron-sized grains of porous amorphous silicates and carbonaceous material surrounded by an important layer of ice, for a 
total dust mass of ~ 0.05M® (in grains up to 1 mm). We discuss evidences that the grains consists of fluffy aggregates. The upper 
limits on the gas atomic lines do not provide unambiguous constraints: only if the PAH abundance is high, the gas mass must be 
;> ' lower than ~ 17M e . 

QO ' Conclusions. Despite the weak constraints on the gas disk, the age of HD 181327 and the properties of the dust disk suggest that 
Q\ , it has passed the stage of gaseous planets formation. The dust reveals a population of icy planetesimals, similar to the primitive 
pT) . Edgeworth-Kuiper Belt, that may be a source for the future delivery of water and volatiles onto forming terrestrial planets. 
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1. Introduction ing the Herschel Space Observatory reveal that as much 
_, ... .as ^30% of nearby Main-Sequence F, G, K stars arc sur- 

The evolution of planets m the Solar System is intimately rounded by cold debris rings ana i g us to the Kuiper Belt 

connected to the existence of a reservoir of planetesimals (lategt He rschel/DUNES Open Time Key Program results, 

i in its outer regions. The depletion of the Kuiper Belt and Eiroa et al in prep) _ Howeverj the cold temperature and the 

g . the mam asteroid belt that occurred during the Late Heavy largc grain sizeg makc - t impossiblc to detcct so lid-state 

Bombardment not only delivered large amounts of water, fcatureg m their spectruni) thus prev enting unambiguous 

volatiles and carbonaceous material onto the inner plan- idcntification of the dust composition. We are thus left with 

ets, but it also reduced the later rate of catastrophic im- model _ depende nt methods to constrain the dust composi- 

pacts onto the Earth, opening the path to the emergence tion from cobr &nd polarimetric measurements, and from 

of life (Morbidelh et al 2000; Greaves & Wyatt 2010; the gpectral E Distribution (hereafter SED). 

Hartogh et al. 2011). Although detecting Kuiper Belt-like . , , r , . . ,. 

, . . iii, , .,1 • Another uncertainty m planet formation theories lies m 

objects around nearby planetary systems will remain an , , , , £ , , ., , / ,. , c 

■' . , . i r .i r iip, ,i • it the lack ol detailed observational constraints for the mu- 

unachievable goal tor the foreseeable future, their colli- , , , , ,. . , . , 

, ° , , „ , r tual dust and gas dissipation tnucscalcs m circumstellar 

sional erosion produces circumstellar disks ot dust rcspon- ,. , . , , ,, , 

., lr I ,-, •• j , , i i 1 ■ disks, or more precisely, on how the gas-to-dust mass ratio 

sible tor a characteristic excess emission detectable at m- , .,, .. , , , . . ,? , , 

r iii -it i ,i i- evolves with time and location m disks at the very early 

irared and sub-milhmctcr wavelengths. Recent studies us- r , , r ^ , £ , , , 

° stages ol planet formation, lhe study ol stellar clusters 

Send offprint requests to: jeremy.lebreton@obs.ujf-grenoble.fr of different ages shows that, statistically, a population of 

* Herschel is an ESA space observatory with science instru- y° un S starS loseS lts maSSlve dust disks m , 0F ^ a few ml1 - 

ments provided by European-led Principal Investigator censor- ll0n y ears on route to the Mam Sequence (Carpenter et al. 

tia and with important participation from NASA. 2005), although, for individual objects, the transition from 

** The Australian Telescope Compact Array is operated by an optically thick to an optically thin disk is expected to oc- 

the Australian Telescope National Facility (ATNF) managed by cur more quickly (a million years or less, Cieza et al. 2007; 

CSIRO. Currie & Sicilia-Aguilar 2011). However not much is known 
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about the characteristic evolutionary timescales for the gas 
in these disks. Observing young stars at different evolution- 
ary stages can help elucidate the extent to which gas and 
dust simultaneously dissipate in disks. This requires obser- 
vations of both the dust and gas components and, at the 
same time, modeling of the continuum and line emission. 

As part of the GASPS (Gas in Protoplanetary 
Systems) Hcrschcl Open Time Key Programme, we ob- 
served HD 181327, a young (~12Myr) F5.5V star located 
at 51.8 pc (Holmberg et al. 2009, Tab. 1) belonging to 
the j3 Pictoris moving group (Zuckcrman & Song 2004; 
Mamajek et al. 2004). HD 181327 was identified as a de- 
bris disk hosting main sequence star through its SED with 
a fractional infrared luminosity Lir/L* = 0.2% (Mannings 
& Barlow 1998). Schneider et al. (2006), from NICMOS 
coronagraphic observations in scattered light, discovered 
a ring-like disk of circumstellar dust located at radius ~ 
89 AU from its star. The ring is inclined by 31.7±1.6° 
from face-on and shows an apparent brightness assymetry 
with respect to the minor axis that is well explained by 
a strong directionally preferential scattering. No photocen- 
ter nor pericenter offset is seen in the ring relative to the 
position of the central star and the disk appears axisym- 
metric. Recent reprocessing of these data (Schneider et al., 
in prep.) confirms the annulus shape of the disk, and sug- 
gests a narrow radial distribution of the dust in the sys- 
tem (FWHM = 24.5 AU). A low surface brightness diffuse 
halo is seen in the NICMOS image and in complementary 
HST/ACS images at a distance of ~ 4" < r < 9" from 
the central star, which may correspond to a population of 
very weakly bound grains originating in the main ring. Mid- 
infrared images at 18.3 //m by Chen et al. (2008) suggest az- 
imuthal asymmetries in the density profile of the ring with 
respect to the apparent minor axis. These 18.3 /mi obser- 
vations also show that little to no emission (< 15%) comes 
from inside the ring resolved in scattered light. Therefore, 
the density profile of the HD 181327 dust disk should be 
well constrained by the scattered light and thermal emis- 
sion images, allowing a detailed study of the dust prop- 
erties. Interestingly, from IRS and MIPS-SED data, Chen 
et al. (2008) also inferred the presence of crystalline ice in 
the debris disk. 

In this paper we perform SED and line fitting of new 
Herschel/PACS observations of HD 181327 (presented in 
Sec. 2), supplemented by data from the literature. We also 
present new ATCA imaging of the disk (Sec. 2). In Sec. 3, 
we derive surface density profiles from the newly reduced 
HST/NICMOS data, which well constrain the ring geom- 
etry. We describe our dust model and present the results 
of the SED fitting in Sec. 3 and 4. We assess the gas con- 
tent in this young debris disk in Sec. 5. Finally we discuss 
the results in Sec. 6, 7 and 8 and present our conclusions in 
Sec. 9. 



2. Observations and Data Reduction 

2.1. Herschel/PACS continuum 

We obtained new far-IR photometry and spectroscopy 
of HD 181327 using the PACS instrument onboard the 
Hcrschcl Space Observatory (Pilbratt et al. 2010; Poglitsch 
et al. 2010). The observations were carried out in two 
modes: 1) 159 sec of point-source chop-nod mode imaging 
at 70 and 160 /im (obsid: 1342183658) and 2) two 276 sec 



Table 1. Star and disk properties 



Parameter 


Value 


Keterence 


Spectral Type 


r 0/r O V 


INordstrom et al. ^zUU4J 


Magnitude V 


/ .U 


ricsscll ^ZUUUJ 


Gravity (log(p)) 




L>nen ex ai. yzuuo ) 


Distance (g£*) 


51.8pc 


Holmberg et al. (2009) 


Age 


12 Myr 


Zuckerman & Song (2004) 
Mamajek et al. (2004) 


Luminosity (L*) 


3.33 l4 a) 


this study 


Mass (M*) 


1.36 Mq 


this study 


Disk PA (b) 


107±2° 


Schneider et al. (2006) 


Disk inclination 


31.7±1.6° 


Schneider et al. (2006) 


im/L* 


~2x 10" 3 


this study 



(a) From integrating the synthetic stellar spectrum, ( ^Position 
Angle, East of North. 



scan map imaging observations at 100 and 160 /im (ob- 
sids: 1342209057-8). The scan maps were executed with the 
medium scan speed (20"/s), at scan angles of 70° and 110°. 
The two cross scan maps at different scan angles were com- 
bined to increase the signal-to-noise by cutting down on 
streaking in the image background along the scan direc- 
tion. 

We reduced the photometry with the Herschcl 
Interactive Processing Environment (HIPE, Ott 2010) ver- 
sion 4.2. The uncertainties in the absolute flux calibration 
for this version are 10% at 70 and 100 /im and 20% at 
160 /im (see the PACS Scan Map Astronomical Observation 
Templates release note 1 ). For the 70" and 160 /an chop- 
node images, we performed aperture photometry using 
apertures of 7" and 11", respectively. For the 100 and 
160 /im scan map observations, we used apertures of radii 
12" and 24", respectively. An annulus for sky background 
estimation was placed from 10" — 20" beyond the aperture 
in each case, and an aperture correction was applied to the 
measured fluxes. Both estimations at 160 /^m are fully com- 
patible within the error bars; the 160 /im scan-map flux is 
used in the rest of the study. Since the source is marginally 
resolved at 70 and 100 /im (see Section 2.1), we calculated 
our own aperture corrections following the procedure de- 
scribed in Poglitsch et al. (2010). The photometry results 
are presented in Table 2. 

Table 2. Photometry results for HD 181327 from chop-nod (1) 
and scan map' 2 ' observations. 



Wavelength 


-*- cont 


Absolute Flux Cal. 


(urn) 


(Jy) 


Uncertainty 




1.827 ±0.0069 


10% 


100 (2) 


1.337 ± 0.0082 


10% 


160 (1) 


0.767 ±0.015 


20% 


160< 2 > 


0.772 ±0.011 


20% 



The ±1<t error bars on the flux measurements include 



statistical uncertainty only. 

HD 181327 is more extended than the Herschel point- 
spread function (PSF) in both the 70 /xm and 100 /im im- 
ages. The 70 /im image of HD 181327 appears in Figure 1, 
along with the 70 /im image of a PSF reference star 

1 available at http://herschel.esac.esa.int/ 

AOTsReleaseStatus . shtml 
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(HD 148387) observed using the same chop-nod mode. The 
upper panel of Figure 2 shows the azimuthally-averaged ra- 
dial brightness profiles for HD 181327 and HD 148387. In 
the case of the 100 (im scan map observation of HD 181327, 
we used a different PSF reference observation, i.e. a scan 
map of a Boo. The 100 /im a Boo radial profile appears 
in Figure 2, bottom panel. At both wavelengths, the az- 
imuthally averaged radial profile of HD 181327 is clearly 
extended beyond that of the PSF reference. The FWHM 
of the HD 181327 70 /Ltm profile is 6'.'44, compared to a 
FWHM of 5'.'28 for HD 148387. At 100 /Ltm, the FWHM for 
HD 181327 is 7'.'58, compared to 6'.'94 for a Boo. 

2.2. ATCA 3.2 mm observations 

The Australia Telescope Compact Array (ATCA) was used 
to observe HD 181327 at 3.2 mm in May and September 
2010. The observations were carried out in continuum dou- 
ble sideband mode, with 2046 x 1 MHz channels per side- 
band with an effective total bandwidth of just over 2 GHz. 
The ATCA has five antenna with 3-mm receivers, and our 
observations were carried out with the compact H75 ar- 
ray configuration with baselines of 31-98 m. The quasars 
PKS 1815-553 and PKS 1253-055 were used for the gain 
and bandpass calibrations respectively, and Uranus was 
used for the absolute flux calibration. The ATCA cali- 
bration uncertainty is estimated to be about 20% in the 
3-mm band. The data were calibrated and reduced with 
the MIRIAD package (Sault et al. 1995). The frequency 
pair for the observations were centered at 93 and 95 GHz, 
and during the reduction process the data was averaged 
over 64 channels to improve the signal-to-noise, shifting 
the frequency pairs to 92.968 and 94.969 GHz. The final 
CLEANed map combined the two sidebands, resulting in 
an image at 93.969 GHz (3192.5 /zm) presented in Fig. 3. 
Our total on-source integration time was 770 minutes. We 
determined an integrated flux of 7.2 x 10 -4 Jy, with an RMS 
of <7 = 1.1 x 10~ 4 Jy/beam. Including the 20% calibration 
uncertainty, this yields an uncertainty of 2.5xlO _4 Jy. The 
disk is marginally resolved. Assuming a Gaussian profile, 
we find a PA of ~ 102°, consistent with that seen on the 
NICMOS image (107 ± 2°, Tab. 1), and along which the 
FWHM is ~ 7.8" (beam: 3.5" x 2.8", PA = 85.4°). The 
center of the Gaussian fit is also shifted by 0.6" toward the 
West with respect to the star position. Structures consistent 




Fig. 1. PACS 70 /xm imaging of (left) HD 181327 and (right) 
a PSF reference star, HD 148387. Both images were acquired in 
chop-nod mode. The field of view of each image is 25" x 25". 
Brightness contours at 5, 10, 25, 50 and 75 times the background 
noise in the HD 181327 image are overlaid on both images. 
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Fig. 2. Radial profiles of HD 181327 at 70 /im (top), compared 
to HD 148387 as a PSF reference, and at 100 /xm (bottom), 
compared to a Boo as a PSF reference. Both plots show the 
azimuthally averaged, normalized radial brightness profile vs. 
radius from the star. Each error bar is the standard deviation 
of the brightness values at that radius. The FWHM values for 
the profiles are given on the plot. 



with a belt at ^100 AU are detected along the semi-minor 
axis with a peak flux of 4.01 x 10~ 4 Jy/beam on the NW 
side. We note a possible brightness asymmetry with respect 
to the SE side. 

2.3. Herschel/PACS Spectroscopy 

We obtained far-infrared spectroscopy using the PACS 
Integral Field Unit which provides a 5 x 5 array of 9'.'4 x 9" 4 
spaxels. The data were taken in chop-nod mode using the 
line scan and range scan modes. The spectra cover eight 
small wavelength ranges (~ 1.8 fj,m wide) targeting partic- 
ular emission lines, including [O I] at 63 /im and [CII] at 
158 fim. The spectral resolution of the data varies between 
R w 1100 and 3400. Additional details on the observations 
appear in Mathews et al. (2010). 

The data were calibrated with the Hcrschcl Interactive 
Processing Environment (HIPE) version 4.2. The 
wavelength-dependent flux and aperture corrections 
provided by the PACS development team were applied, 
giving an uncertainty in the absolute flux calibration 
of 30% (see the PACS Spectroscopy Performance and 
Calibration Document 1 ). Significant emission was detected 
only in the central spaxel. No emission lines are seen, only 
continuum flux. The central spaxel spectrum centered on 
the [O I] 63 /Ltm line appears in Figure 4. 

The continua were binned at half the instrumental 
FWHM and fitted with 1st degree polynomials, using the 
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RA (arcsec) 

Fig. 3. ATCA 93.969 GHz (3.192 mm) map of HD 181327 with 
contours at -2 (dashed red line), 2 and 3a (solid blue lines), 
where a is the image RMS given by l.lxlCP 4 Jy/beam. The 
dashed blue area represents the beam. The Gaussian fit is over- 
plot with contours at 2 and 3a (solid green lines). The red 
and green crosses denote respectively the star position and the 
Gaussian fit center. 



fluxes and errors produced by our HIPE reduction script. 
The HIPE uncertainties appeared to be slightly too small, 
however. New flux random errors were calculated from the 
RMS of the data minus the continuum fits, over wavelength 
ranges spanning 6 to 14 times the instrumental FWHM, 
centered on the expected line rest wavelengths. The contin- 
uum fluxes at the rest wavelengths of the expected emission 
lines appear in Table 3. To determine the upper limit on 
the integrated line flux (-Pint) for each available emission 
line, we integrated the data over a small wavelength range, 
either ±1 or ±1.5 resolution elements, centered on the rest 
wavelength of the expected emission line. The flux errors 
were propagated to determine the la upper limit on each 
integrated line flux. The 3er integrated line flux upper limits 
appear in Table 3 as well. 



Table 3. Spectroscopy Results for HD 181327 



~3 



1.5 




Line 


Wavelength 


cont 


1 int 




(H 


(Jy) 


(xl(T 18 W/m 2 ) 


LoiJ 


63.1852 


1.731 ±0.236 


< 9.82 


CO 36-35 


72.8429 


1.491 ±0.185 


< 18.41 


0-H2O 


78.7414 


1.431 ±0.187 


< 10.51 


CO 29-28 


90.1630 


1.264 ± 0.222 


< 8.90 


[01] 


145.5350 


0.804 ±0.107 


< 8.48 


[CII] 


157.6790 


0.767 ±0.175 


< 7.02 


0-H2O 


179.5265 


0.388 ± 0.201 


< 6.87 


DCO+ 22-21 


189.5700 


0.581 ± 0.378 


< 16.71 



^ The ±lcr error bars on the flux measurements include sta- 
tistical uncertainty only. The absolute calibration uncertainty 
amounts to 30%. * 3a upper limit 



63.0 63.1 63.2 63.3 63.4 
Wavelength (urn) 

Fig. 4. PACS spectrum of HD181327, centered on the wave- 
length of the [O I] 63 fj,m emission line. The linear fit to the 
continuum flux is overplotted with the red line. The [O I] 63 /J,m 
emission line was not detected. 



3. Model Setup 

The dust and gas models developed in Sees. 4 and 5 rely 
on some assumptions on the spatial distribution of the dust 
and the gas, as well as on the dust properties. These are de- 
tailed below in Sees. 3.1 and 3.4. Additional data collected 
from the literature and stellar properties are summarized 
in Sec. 3.2. 



3.1. New Surface Density Profile for the HD 181327 Ring 

Our spatially resolved imaging constrains the location of 
the dust. We advantageously impose this constraint in order 
to examine the dust properties in details and to get an up- 
per limit on the gas mass from the Herschel non-detections 
of gaseous emission lines. The disk around HD 181327 was 
spatially resolved in scattered light with HST/NICMOS 
1.1 //in imaging by Schneider et al. (2006), which alleviates 
many uncertainties on the disk structure. The NICMOS im- 
ages revealed a ring-like disk with a peak density position 
ro = 1.71" (~ 88.6 AU, according to the new distance to the 
star, 51.8 pc, by Holmberg et al. 2009), with an azimuthally 
medianed photometric FWHM of - 0.538" (~ 28 AU). 

Recently, Schneider et al. (in prep.) improved the data 
reduction of the original HST/NICMOS observations of 
HD 181327 thanks to PSF subtraction template images bet- 
ter matched to the stellar PSF structure in the HD 181327 
images. This results in an overall improvement of the signal- 
to-noise, especially inside the ring, and leads to a slight 
sharpening of the belt. The new profiles, displayed in the 
top panel of Fig. 5, show a peak at a distance of 88.6 AU and 
a surface brightness fainter by ~ 21% (0.76 mJy. arcsec ~ 2 at 
the peak position) compared to the original Schneider et al. 
(2006) study. 

We derive the surface density profile S(r) by invert- 
ing the new scattered light surface brightness profile. We 
employ the methodology described in Augereau & Beust 
(2006), assuming a disk inclination of i = 31.7° from face- 
on. The inversion procedure yields the product <7 sca £(r), 
where C7 sca is the mean scattering cross section of the grains 
at 1.1 fim. This is equivalent to the product of the ver- 
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Fig. 5. New HST/NICMOS 1.1 /im surface brightness profiles 
(top, nominal and ±1ct profiles), and surface density profiles 
resulting from direct inversion of the brightness profiles (bot- 
tom), assuming two possible values for the asymmetry param- 
eter. The color-coding is the same for both panels. The striped 
areas identify regions dominated by PSF-substraction residuals 
where the density has been extrapolated. 



ported by the mid-IR results of Chen et al. (2008), who 
show that the inner part of the disk is largely cleared to- 
ward the center, containing little to no dust. 

Outside of the belt peak position, the improved data 
reduction results in a profile that falls off as ~ r -4.7±o.3 
from ~92 to 162 AU, a steeper fall-off than the Schneider 
et al. (2006) original profiles (~ r~ 3 ). The dust is thus 
more tightly confined to the parent-belt than previously 
estimated, and the profile is also steeper than theoretical 
expectations. The exact value of the outer slope can only 
slightly impact the modelling results, as most of the emis- 
sion arises from the inner part of the belt where the temper- 
ature is larger, except possibly for the longer-wavelengths. 

Eventually, the differences between the several radial 
profiles in Fig. 5 essentially translate into a scaling factor 
that impacts the total disk mass estimate. We estimate that 
the uncertainty on the asymmetry parameter adds an ad- 
ditional ^4% to the uncertainty on our final mass values. 
The density profile obtained for |<?hg |=0.30 is used in the 
rest of this study. We stress that, in order to focus on the 
other disk properties, we make the assumption that this 
profile is valid for any dust grain constrained by the SED, 
i.e. for grain sizes ranging from the micron to the millime- 
tre scale. This assumption assumes no radial segregation in 
dust sizes and it remains reasonable as long as the flux at 
any wavelength arises predominantly from regions close to 
the density peak position (see Sec. 6.2). 



tical optical depth r±(r), multiplied by the mean albedo 
uj at that wavelength (<x 8ca £(r) = u)T_\_(r), see Eqs. 2-3 in 
Schneider et al. 2006). The inversion solution depends on 
both the anisotropic scattering phase function, and on the 
disk vertical profile. Nevertheless, Schneider et al. (2006) 
showed that the density profile is largely independent of 
the shape for the vertical dust distribution in the range 
of acceptable vertical scale heights Hq. We here assume a 
Gaussian vertical profile with a scale height Ho = 6AU 
at the observed peak radius of the ring r = 88.6 AU 
(Ho/ro ~ 0.068, consistent with the theoretical lower limit 
of 0.04 from Thebault 2009). The observed phase func- 
tion of the disk at 1.1 (im is well-represented by a Henyey- 
Greenstein scattering phase function with asymmetry pa- 
rameter |<7hgI= 0.30 ± 0.03. This value was discussed by 
Chen et al. (2008) in light of asymmetries identified in T- 
RcCS thermal emission images in the Q a -band. Chen et al. 
(2008) suggest that the asymmetries may have led to an 
overestimate of <?hg|, which may be as small as 0.16±0.04. 
We thus perform the inversion for both |<?hg| = 0.16 and 
\9hg\= 0.30. 

The results are displayed in the bottom panel of Fig. 5, 
where we also show the impact of the surface brightness un- 
certainties on the inferred surface density profile. We find 
an actual surface density peak at 89.5 AU from the star. 
A higher </hg leads to a globally higher surface density 
profile. The surface brightness measured interior to ~ 1.5" 
(~ 78 AU) is not reliable because of PSF subtraction resid- 
uals. We fit the profiles from 79 to 83 AU obtained for both 
|.9hg| values by a very steep power-law function of the form 
~ r +10±1 . We extrapolate this profile interior to 78 AU as 
shown in Fig. 5 (striped area). This extrapolation is sup- 



3.2. Spectral Energy Distribution 

The Herschel and ATCA data presented in this paper are 
complemented by IR and sub-mm data collected from the 
literature. These are listed in Table 4. The shortest wave- 
length flux we use is the 12 /im IRAS photometry for which 
the emission shows no excess from the stellar photosphere 
level. This observation is documented in Schneider et al. 
(2006) together with IR-excess from the IRAS Faint Source 
Catalogue and obtained with the Spitzer/MIPS instrument. 
These are complemented by ground-based mid-IR measure- 
ments performed with the T-RcCS imager on Gemini-S by 
Chen et al. (2008), and 13 points extracted from the 5- 
35 ^m Spitzer/IRS spectrum that we reduced using the c2d 
Legacy team pipeline (Lahuis et al. 2006). 

Additional MIPS observations were obtained in SED- 
mode by Chen et al. (2006) between 54.5 and 95.5 /im, 
well-constraining the shape of the disk SED around its max- 
imum emission wavelength. Far-IR excess was reported by 
Moor et al. (2006) using the ISO database. The shape of the 
SED in the far-IR domain (100-200 /im) is well constrained 
by our new Herschel measurements. 

Because the uncertainties on the dust optical constants 
at millimetre wavelengths arc high, we chose not to use 
the ATCA flux to constrain the dust model (see Sec. 6.2 
for a discussion on the compatibility of the model to these 
observations). Therefore, the sub-mm regime is only rep- 
resented by a measurement with the LABOCA bolometer 
array at the 12-m telescope APEX (Nilsson et al. 2009, 
A ~ 870 /im), which thus represents a key piece of informa- 
tion. These archival data, together with our new observa- 
tions, constitute the 49 measurements that are displayed in 
Fig. 6. 
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Table 4. HD 181327 photometric observations from the litera- 
ture. The Herschel continuum values used for the dust modeling 
are in Tables 2 and 3. 



Wavelength 


Flux 


Uncertainty 


Instrument 


(cm) 


(Jy) 


(Jy) 


& Reference 


12 


0.114 


0.011 


IRAS [1J 


15.0 


0.094 


0.004 


O " J- /TT> o roi 

Spitzer/IRS [2J 


16.3 


0.097 


0.005 


Spitzer/IRS [2] 


17.4 


0.107 


0.005 


O 'J- /TT> o roi 

Spitzer/IRS [2J 


18.3 


0.114 


0.003 


T-ReCS [3J 


18.7 


0.125 


0.007 


O " J- /TT> O Pol 

Spitzer/IRS [2J 


19.8 


0.144 


0.009 


Spitzer/IRS [2J 


20.4 


0.149 


0.008 


Spitzer/IRS [2J 


21.1 


0.157 


0.020 


Spitzer/IRS [2] 


23.1 


0.224 


0.007 


Spitzer/IRS [2] 


23.7 


0.223 


0.009 


Spitzer/MIPS [1] 


25 


0.306 


0.024 


IRAS [1] 


25.7 


0.303 


0.014 


Spitzer/IRS [2] 


28.0 


0.384 


0.014 


Spitzer/IRS [2] 


30.6 


0.488 


0.026 


Spitzer/IRS |2J 


33.1 


0.623 


0.025 


Spitzer/IRS [2] 


35.5 


0.718 


0.034 


Spitzer/IRS [2] 


54.5 


1.88 


0.09 


MIPS-SED [4] 


58.0 


1.83 


0.07 


MIPS-SED [4] 


60 


1.73 


0.17 


ISO [5] 


61.4 


1.98 


0.08 


MIPS-SED [4] 


64.9 


2.01 


0.08 


MIPS-SED [4] 


68.4 


1.98 


0.08 


MIPS-SED [4] 


71.4 


1.73 


0.12 


Spitzer/MIPS [1] 


71.6 


1.99 


0.09 


MIPS-SED [4] 


75.1 


1.93 


0.10 


MIPS-SED |4J 


78.4 


1.81 


0.10 


MIPS-SED [4] 


81.8 


1.72 


0.11 


tv TTn t~i i — * r a 1 

MIPS-SED [4] 


oo.o 


1 R9 


n 1 1 


MTP^ SFD T4l 
iviiir o-O-Cj-L' [4:1 


88.8 


1.47 


0.11 


MIPS-SED [4] 


90 


1.41 


0.14 


ISO [5] 


92.0 


1.46 


0.13 


MIPS-SED [4] 


95.5 


1.37 


0.17 


MIPS-SED [4] 


100 


1.7 


0.2 


IRAS [1] 


155.9 


0.77 


0.09 


Spitzer/MIPS [1] 


170 


0.736 


0.192 


ISO [5] 


870 


0.052 


0.006 


LABOCA [6] 


3190 


7.2xl0~ 4 


2.5xl0~ 4 


ATCA [7] 



Notes and References - All uncertainties are la and include 
absolute calibration error. [1] Schneider et al. (2006), [2] This 
study, data reduced with the Spitzer/c2d team pipeline (Lahuis 
et al. 2006), [3] Gemini, Chen et al. (2008), [4] Spitzer, Chen 
et al. (2006), [5] Moor et al. (2006), [6] APEX, Nilsson et al. 
(2009), [7] This study, Sec. 2.2. 



3,3. Star Properties and dik geometry 

We derive excesses from a comparison of the SED with a 
synthetic stellar spectrum extracted from the "NextGcn" 
model atmosphere grid (Hauschildt et al. 1999) and scaled 
to match the visible magnitude V = 7.0 mag (Bcssell 2000), 
assuming log(g) = 4.5, T c ft = 6600K, which leads to a stel- 
lar luminosity of 3.33 L© and a mass of 1.36 M Q . This stel- 
lar model shows good agreement with the Spitzer/IRS spec- 
trum at wavelengths shorter than 10-15 /im, longward of 
which the excess starts dominating over the stellar pho- 
tosphcric emission. The star properties and disk geometry 
adopted in this study are summarized in Table 1. 



10.000 




1 10 100 1000 

Wavelength [|jm] 



Fig. 6. Observed Spectral Energy Distribution of HD 181327 
(see Tables 2, 3 and 4) 



3.4. Dust Properties: Size Distribution and Composition 

For both the dust and gas modeling, we make the assump- 
tion that the grain population has a single differential size 
distribution valid anywhere in the ring. This is a central 
hypothesis that necessarily limits the uniqueness of the so- 
lution. We discuss its impact in Sec. 6. The differential size 
distribution follows a classical power law dn(a) oc a K da, 
from a minimum size a m i n to a maximum size a max - The 
minimum grain size a m i n and the power law exponent k 
are free parameters that are constrained by the SED fit- 
ting (Sec. 4). The maximum grain size a max is a fixed pa- 
rameter, sufficiently large so that it does not affect the 
light emission properties at the wavelengths considered (we 
chose a max ~ 8 mm). For standard k values (close to -3.5 or 
smaller), the mass of solid material in a debris disk can be 
"arbitrarily" large depending on the assumed higher cut- 
off value of the size distribution (that can extend up to 
planetcsimal sizes). In this article, we define the dust mass 
-Mdust as the total mass of material contained in grains 
smaller than 1 mm. This definition has been followed within 
the GASPS consortium because it allows direct comparison 
between different studies independently from the model as- 
sumptions. 

We adopt a complex grain model that bears some sim- 
ilarities to those succesfully employed by several authors 
for fitting SEDs of young debris disks (Li & Grecnbcrg 
1998; Augereau et al. 1999; Li & Lunine 2003a,b), and 
for fitting scattered light disk colors (Kohler et al. 2008). 
We consider icy porous aggregates, with silicate cores 
and organic refractory mantles (or simplified versions of 
this model; see Sec. 4), to mimic fluffy particles made of 
essentially unalterred interstellar material (the so-called 
"cold-coagulation" dust model in Li & Lunine 2003a, b). 
Assuming cosmic abundances, the volume of the carbona- 
ceous mantle Vq is estimated to be ~ 2.12 times that of the 
silicate core Vsi in the interstellar medium, a value similar 
to the one derived from the mass spectra of comet Hallcy 
in the Solar System (V c /V si - 2, Greenberg 1998). 

To calculate the optical properties, we use a num- 
ber of representative materials with well established op- 
tical indexes. These are either derived from laboratory 
measurements or from observational and analytical argu- 
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Table 5. Parameter space explored with GRaTer, corresponding to a total of about 1,150,000 models. The volume fractions vsi, 
dc, «ic C and the porosity V are denned in Sec. 3.4 (Eq. 1) 





Description 


Parameter 


Explored range 


Values 


Distribution 


Sec. 3.4 


minimum grain size 
size power law index 
maximum grain size 
Disk mass up to 1 mm 


a m in [/xm] 

amax [mm] 
M dust [M e ] 


0.01... 100 
-6.0. ..-2.5 
7.8 

> 


77 

50 
fixed 
fitted 


log 
linear 


1 AND 2-MATERIAL MODELS 


Sec. 4.1.1 
Sec. 4.1.2 
Sec. 4.1.2 
Sec. 4.1.3 


Si+ Amorphous C 
Si+amorphous ice 
Si+crystalline ice 
Porous Si 


ACAR volume fraction vc 
ice volume fraction Vi ce 
ice volume fraction Vi ce 
Porosity V 


0.00. . . 1.00 
0.00 ...0.90 
0.00 ...0.90 
0.00 . . . 0.99375 (a) 


21 
19 
19 
24 


linear 
linear 
linear 
linear 


3 AND 4- MATERIAL MODELS 


Sec. 4.2 


Si 

+amorphous C 
+ amorphous ice 
+vacuum 


astrosi volume fraction Dsi 
am. carbon volume fraction vc 
ice volume fraction Vi ce 
Porosity V 


l/3x(l-« ice ) 
2/3x(l-« ice ) 
0.00. ..0.90 (b) 
0.00. ..0.95 (c) 


12 
18 


linear 
linear 



Notes - (a) 19 linearly distributed points from 0.00 to 0.95 + additional high porosity points : 0.96875, 0.975, 0.9875, 0.99375. 
' b ' 10 linearly distributed points from 0.00 to 0.90 + additional points at 0.65 and 0.75. ' c ' 10 linearly distributed points from 0.00 
to 0.90 + additional points at 0.45, 0.55, 0.625, 0.65, 0.675, 0.75, 0.85, 0.95. 
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10.0 



1.0 



0.1 



\ .-. 

\ 

:7::r,,.,A^ 

Compact Si (p=3.50g/cm 3 , x 2 =1 0.3) 
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Fig. 7. Mass opacity of spherical grains computed with the Mie 
theory and Bruggeman mixing rule for several prototypical com- 
positions detailed in Tab. 6 and their associated best-fit size dis- 
tributions. 



close to ellipsoids such that the porosity can effectively be 
seen as the volume fraction of vacuum (Shen et al. 2009, 
and references therein). In fact, Shen et al. (2009) show 
that only the polarization and the phase functions are sig- 
nificantly sensitive to grain shapes. Wc define the composi- 
tion of a grain through the volume fractions of its different 
components : 

Vsi Vq Vi ce Kac / 1 \ 

V Si = Tr , V C = — , V icc = — , r = — — (1) 



v m 



v m 



v m 



v Vi 



where Vsi, Vc, Vi cc and V vac are the total volume of astro- 
nomical silicates, carbonaceous material, ice and vacuum 
in the grain, respectively, and V ma t = Vsi + Vc + Mco is 
the total volume of solid material. V denotes the poros- 
ity. Should one of the materials within the grains reach its 
sublimation temperature, it is replaced by vacuum and the 
optical properties are adjusted accordingly in the model. In 
most cases however, the disk is too cold for any material 
to sublimate. Fig. 3.4 shows the mass opacity (absorption 
cross section per unit mass) calculated using Eq. 4 of Draine 
(2006) for several prototypical dust compositions, with the 
size distributions given in Tab. 6. 



ments. The materials are: (1) the astronomical silicates (as- 
troSi) from Draine (2003), a theoretical material consist- 
ing of amorphous silicates polluted by carbonaceous ma- 
terials and other metals (p = 3.5 g.crn" 3 ); (2) the ACAR 
sample of amorphous carbon from Zubko et al. (1996) 
as organic refractories (p = 1.95 g.cnr 3 ); (3) the H2O- 
dominated ice from Li & Greenberg (1998) (amorphous 
ice, p = 1.2 g.cm -3 ); and (4) the crystalline water ice from 
Warren (1984) (p = 0.92 g.cm" 3 ). 

The optical indexes of an aggregate are obtained using a 
standard effective medium theory (EMT), the Bruggeman 
mixing rule (Bohren & Huffman 1983), and following the 
methodology of Li & Greenberg (1998) and Augereau et al. 
(1999). The optical properties are computed with Mie the- 
ory valid for hard spheres. The EMT-Mic combination has 
been shown to be a good approximation for the absorption 
and scattering efficiencies, provided that grain shapes are 



4. Modeling the Dust Continuum Emission 

We use the GRaTer code (Augereau et al. 1999) to com- 
pute a grid of models with a goal of reproducing the 
HD 181327 Spectral Energy Distribution. GRaTer is a 
SED, image and interferometric data fitter, specially de- 
signed to efficiently model optically thin disks with para- 
metric grain size and radial distributions, or distributions 
from dynamical simulations. The code accounts for both 
the scattered light and continuum emission of dust grains 
in thermal equilibrium with a star. 

In the present case, the spatial dust distribution is con- 
strained by the resolved HST images (Sec. 3.1, Fig. 5). The 
surface density profile derived from inverting the observed 
surface brightness was imposed in the model. Our study of 
the disk continuum emission focuses on the properties of 
the dust particles, namely their size distribution and com- 
position. To that end, we computed large grids of models 
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Fig. 9. Bayesian probabilities for the models in Table 6: minimum grain size, size distribution power law index, and total dust 
mass. 



1.00 


/%! 


A 

\ \ 

\ \ 


0.10 


■ Compact Si (y 2 =1 0.3) 

Si +\£ (x 2 =775) , 

Si + amorph. ice (x =8.90) 

Si + cryst. ice (v =10.4) 

Porous » (x =6.05) 


X \ 

V V 

\ H*. 


0.01 






10 100 1000 
Wavelength [um] 


1.00 


\ // 
■ \ if 1 




0.10 




\ 

\ V _ 
\ V 
\ V 

\ ] 




Si+C\Porosity (x 2 =6.05) \ \ 
--- Si+C+ice (x 2 =4.48) \ X 
Si+C+ice+Porosity (x =1 .64) \ \ 


0.01 




i \\ 



10 100 1000 

Wavelength [um] 

Fig. 8. Best fit models (Table 6) for several 1- or 2-material 
mixtures (top) and 3- to 4-material mixtures (bottom). Red 
crosses: photometric data (Fig. 6). Solid grey line: synthetic stel- 
lar spectrum. 

reflecting various possible dust properties, starting from the 
simplest composition (astronomical silicates) and progres- 
sively increasing the model complexity until obtaining a 
good match to the SED. For each model, the dust mass 
•Mdust and the quality of the fitting are obtained from lin- 
ear least-square fit between the observed and the synthetic 
SEDs. The x 2 listed in this section are obtained using the 
48 data points presented in Tables 2, 3 and 4, and with the 
numbers of degrees of freedom listed in Table 6. We sub- 
sequently use a statistical (Bayesian) inference method to 



derive the probability distribution of each parameter (see 
Appendix A). 

In total we consider ~ 300 compositions, with 77 x 50 
grain size distributions for each, leading to ~ 1,150,000 
models. All the parameters of the grain models are sum- 
marized in Table 5. 

4.1. First Steps Toward a Complex Dust Model 

We start by fitting the HD 181327 SED with a single- or 
two-material grain model, which show the need for a more 
complex dust model. We consider several mixtures of the 
materials presented in Sec. 3.4 to assess the broad proper- 
ties of the grains. The model then includes 3 or 4 free pa- 
rameters: (1) the minimum grain size a m ; n ; (2) the slope of 
the size distribution /c; (3) the total mass of the disk Afd us t; 
and in the case of 2-component grains: (4) the volume frac- 
tion of either carbon (vc) or ice (Vi ce ), or the porosity (V). 
The detailed results of the investigation are summarized in 
Table 6, the best-fitting SEDs are displayed in Fig. 8 (top 
panel) and probability curves are shown in Fig. 9. 

4.1.1. Silicates or Organic Refractories? 

The most basic grain model to consider is the case of com- 
pact spheres made of pure astronomical silicates (astroSi), 
a configuration that is often used by default in many model- 
ing approaches. We let a m - ln , k and Mdust vary and we com- 
pute a grid of x 2 on which we apply a Bayesian analysis. 
We find that pure astronomical silicate grains would have 
Omin = 0.99 ± 0.13 /im and k = —3.43 ± 0.04, leading to a 
total dust mass Md us t = 0.35 ± 0.03M® for grains smaller 
than 1mm (Tab. 6). These values are consistent with the 
ones typically found for debris disks, but they actually 
lead to very poor fits, as illustrated in Figure 8 (top panel) 
where the model with the smallest x 2 {Xi = 10.3) is dis- 
played. Although the model is in fairly good agreement 
with the data over the full spectral range and agrees well 
with the sub-mm measurement, it fails to reproduce the 
shape of the emission peak at ~ 50 — 100 (im and the 
mid-IR fluxes, shifting the dominant part of the disk to 
colder temperatures. It overestimates by 30-40% the excess 
at far-IR wavelengths (A ~ 100 — 160 /mi), in particular the 
Herschel/PACS measurements. 

We take a further step and consider mixtures of amor- 
phous silicates and carbonaceous materials. We ran a 
Bayesian analysis and find that the actual best composi- 
tion corresponds to pure carbon. Qualitatively, the effect 
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Table 6. Best-fitting (Best, i.e. lowest Xr) an d most probable (Proba, i.e. highest probability) parameters assuming different grain 
compositions (see text for the definition of the parameters). The most probable parameters are given with 3-cr uncertainties. The 
volume fractions marked with the I symbol are imposed knowing the other free volume fractions. 



Composition 




VSi 


vc 


Vice 


V 


a m i n (/xm) 


K 


M du8t (M e ) 


xi 


1-MATERIAL (dof = 45) 


Si 


Proba : 
Best : 


1.00 
1.00 








0.99±0.13 
1.00 


-3.43±0.04 
-3.42 


0.35±0.03 
0.35 


10.3 


C 


Proba : 
Best : 




1.00 
1.00 






1.62±0.03 
1.62 


-3.65±0.07 
-3.64 


0.15±0.03 
0.15 


7.75 


2-MATERIAL (dof = 44) 


Si+C w 


Proba : 


0.0071 


0.99±0.04 






1.62±0.03 


-3.65±0.07 


0.15±0.03 




Si+am.ice 


Proba : 
Best : 


0.591 
0.60 




0.41±0.05 
0.40 




0.15±0.13 
0.16 


-3.36±0.07 
-3.35 


0.26±0.05 
0.26 


8.90 


Si+cr.ice 


Proba : 
Best : 


0.67f 
0.65 




0.33±0.12 
0.35 




0.29±0.24 
0.26 


-3.35±0.02 
-3.35 


0.27±0.04 
0.27 


10.4 


Si+P 


Proba : 
Best : 


0.941 
0.025 






0.057±0.062 
0.975 


0.98±0.20 
0.010 


-3.42±0.05 
-3.28 


0.35±0.04 
0.039 


11.9 


3-MATERIAL (dof = 44) 


Si+C+V 


Proba : 
Best : 


0.33+ 
0.331 


0.67+ 
0.67f 




0.95±6xl0" b 
0.95 


9.45±3.20 
10.0 


-3.13±0.09 
-3.14 


0.042±0.001 
0.042 


6.05 


Si+C+ice 


Proba : 
Best : 


0.07f 
0.07f 


0.13f 
0.13+! 


0.700±0.002 
0.80 




0.78±0.04 
0.55 


-3.50±0.03 
-3.42 


O.lliO.01 
0.11 


4.48 


4-MATERIAL (dof = 43) 


Si+C+ice+P 


Proba : 
Best : 


O.llf 
0.12f 


0.22f 
0.23+! 


0.67±0.07 
0.65 


0.63±0.21 
0.65 


0.81±0.31 
0.89 


-3.41±0.09 
-3.42 


0.051±0.016 
0.048 


1.64 



Notes - "dof" means degrees of freedom (see Appendix A). The best fit for the Si+C mixture is obtained for silicate-free 
grains. The best-fitting parameters for that solution are given in the 1-material part of the table. 



on the SED is not significant (see Fig. 8, top panel), even 
if it leads to a better minimum chi-square (x 2 = 7.75). 
This model gets closer to the measured fluxes at the peak 
emission (A ~ 50 — 100 /im), but does not provide an over- 
all good fit. The addition of carbonaceous material slightly 
shifts the best fit parameters of the grain size distribution 
(Tab. 6), resulting in slightly larger a m ; n values and slightly 
steeper size distributions 



4.1.2. Icy Grains? 

Chen et al. (2008) suggested that HD 181327 may contain 
a significant fraction of water in the form of crystalline ice, 
producing a broad peak in the combined Spitzer/IRS (5- 
35 /an) and MIPS-SED (60-77 /im) spectrum. To assess the 
presence of water ice in the debris disk of HD 181327, we 
model the disk with grains made of a mixture of astronom- 
ical silicates and water ice, either in the form of amorphous 
ice or crystalline ice. 

The results of the statistical analysis are documented 
in Table 6. The quality of the fit is only slightly improved 
compared to the pure silicate model: the best models have 

= 8.90 for amorphous ice, Xi = 10-4 for crystalline ice. 
In fact, the addition of water ice, crystalline or amorphous, 
does not modify the general shape of the synthetic SED (see 
Fig. 8, top panel). The smallest x 2 calculated with amor- 
phous ice is slightly smaller than the one of crystalline ice, 
but in both cases the values are too high to accept these 
models. 



4.1.3. Porous Grains? 

As a final step in our preliminary study of the dust com- 
position, we consider a grid of models consisting of sili- 
cate grains with different porosities, from compact silicates 
to extremely fluffy silicate grains (P max = 99.375%). The 
marginal probability distribution for the porosity V re- 
veals a bimodal probability distribution: the most probable 
models arc the ones with either no porosity (discussed in 
Sec. 4.1.1), or very high porosity (V = 97.5%). Both solu- 
tions are presented in Tab. 6. A visual inspection of Figure 8 
resolves the ambiguity. The compact grain model better 
reproduces the minimum and maximum wavelength ends 
of the SED, and in particular fits the sub-mm point, but 
it fails at reproducing the intermediate wavelength range 
around the emission peak. On the other hand, the highly 
porous grain model is orders of magnitude below the ob- 
served value at 870 /im, but it is within the error bars in 
the other regions of the SED, especially at Herschel wave- 
lengths. In sum, none of these models lead to overall good 
fits to the whole SED. 



4.2, Construction of an Advanced Grain Model 

From what precedes, we conclude that none of the 2- 
component mixtures considered are sufficient to provide 
a satisfactory fit to the SED. This is essentially a conse- 
quence of the disk surface density being known a priori in 
this model; it cannot be adjusted to compensate for the 
emission properties of the dust grains. We note, however, 
that the addition of one component to silicates, should it 
be carbonaceous material, amorphous ice, or porosity, pro- 
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Fig. 10. Bayesian probabilities for the grain composition in the case of 4-material dust mixtures. 



vides some slight improvements. This suggests that their 
combined effect could produce a more realistic model. 

We evaluate the improvement provided by a 3- 
component mixture on the fitting of the SED. We showed in 
Sec. 4.1.1 that it is not clear how to disentangle the respec- 
tive contribution of silicates and organic materials by look- 
ing at the SED. Instead of including an additional free pa- 
rameter, we consider grains (approximately) in fixed Si/C 
ratio described in Sec. 3.4 (vc/vsi — 2), and we restrict the 
discussion to the ice fraction and the porosity of the grains. 
We thus assume that grains are made of silicate and car- 
bonaceous material, which are either icy (wi C0 parameter) or 
porous (V parameter). The results are reported in Tab. 6. 
For porous ice-free grains, the best model is found for a 
porosity V — 95% and leads to some improvement to the fit 
with = 6.05. If we now consider only compact icy grains, 
we can achieve a very reasonable x 2 - = 4.48 with compari- 
son to previous models, assuming an amorphous ice fraction 
Vice = 0.7. Overall, both the addition of either water ice or 
vacuum to silicate + carbonaceous grains do contribute to 
improving the fit to the SED of HD 181327. This can can 
be seen in Figure 8 (bottom panel). Nevertheless, the fits 
are still unsatisfactory (x 2 much larger than 1), and justify 
an even more advanced dust model. 

We next consider a 4-component grain model inspired 
by the "cold coagulation dust model" by Li & Lunine 
(2003a, b), which has been successfully employed for model- 
ing debris disks (Li & Greenberg 1998; Augereau ct al. 1999; 
Li & Lunine 2003a,b; Kohler et al. 2008). We keep the Si/C 
volume ratio fixed (usi:"Uc=l:2) and let both the ice fraction 
and the porosity vary, in addition to the other usual free 
parameters in this study, namely the minimum grain size 
cfmim the slope of the size distribution k, and the total mass 
of the disk Md us t- This leads to a five- free-parameter disk 
model to fit the 48 data points in the SED. We compute 
three grids of models corresponding to the three surface 
density profiles shown in Fig. 5 for |<7hg| = 0.30 (i.e. the 
nominal and ±lcr profiles). The Bayesian analysis is then 
applied to the 3 grids separately and the probabilities for 
each point of the parameter space are calculated as the sum 
of the probabilities for the 3 profiles. The detailed results 
are documented in Tab. 6. 

Figure 10 presents the results of the Bayesian analysis 
on the composition of the dust grains, with a 2D- map show- 
ing the probability as a function of V and Vi ce for the grid 
models, after projection of the probabilities over the other 
parameters. It clearly illustrates the need for a composition 
within a limited range of possible values. In particular, the 
fraction of H20-dominated ices is well defined: the proba- 



bility peak is obtained for an ice fraction Wi CC = 0.67 ±0.07. 
The porosity is not so strongly constrained but we find a 
peak of probability for V = 0.63 ± 0.21. The best fit model 
in terms of x 2 is found for grains made of 23% carbona- 
ceous materials, 12% silicates, 65% ice, and with a porosity 
V = 65%. As can be seen in the bottom panel of Figure 8, 
for such grains, the synthetic SED very reliably describes 
the observed SED from the mid-IR to the sub-mm, with 
X 2 = 1.64, a value significantly better than those obtained 
for simpler grain models. The grain size distribution de- 
rived from the statistical analysis has a m ; n = 0.81 ± 0.31 
/jm and n = —3.41 ± 0.09, a power law index very con- 
sistent with that predicted for theoretical collisional cas- 
cades (k = —3.5). The mass in grains smaller than 1 mm, 
-Mdust = 0.051 ± 0.016 Mq, may look surprisingly small 
at first glance compared to the other models, but it re- 
sults from the significant fractions of ice and vacuum in the 
grains, implying larger mass opacities (Fig. 3.4). The tem- 
perature of the grains at the surface density peak position 
(89.5 AU) ranges from 86 K for the smallest grains to 40 K 
for ~ 8 mm-sized grains. 

To summarize, our approach, based on a fixed surface 
density profile constrained by resolved imagery of the disk, 
clearly highlights the need to consider the presence of both 
ice and porosity to find an overall good fit to the SED. 
Icy grains are naturally expected at such a large distance 
from a solar-type star and the fluffincss of the grains is 
likely a by-product of their formation process. We applied a 
similar cold-coagulation dust model as Li & Lunine (2003a) 
did fitting the very young HD 141569A debris disk, except 
that we used no a priori assumption on the ice fraction, 
but with arguing the grain composition must reflect the 
primordial abundances of condensable elements, they chose 
ti ico ss 30% for cold regions and they find a best-fitting 
porosity V « 0.8. Hence we propose a dust model that 
is especially ice-rich with a degree of porosity that is not 
unexpected. Further comparisons with the literature and 
implications on the nature of the dust are discussed in more 
details in Sec. 6. 



5. Limits on the Gas Content 

At 12 Myr-old, HD 181327 offers a chance to assess the 
amount of gas in a young debris disk, at a time when, ac- 
cording to planet formation theories, the accretion of gas 
by giant planets has already occurred while the formation 
of rocky planets might be ongoing. 

We use the upper limits on the oxygen and ionized 
carbon fine-structure lines ([O I] at 63 and 145 /im, and 
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[C II] at 158 fj,m, see Tab. 3) to constrain the amount of 
gas present in the disk. We solve the photochemistry and 
energy balance in the disk using the photochemical code 
ProDiMo 2 . The abundance of 71 species is computed in 
steady-state, under the influence of adsorption, desorption 
(thermal, photo and cosmic-ray induced), and molecular 
hydrogen formation on grain surfaces. 

The gas density profile is supposed to match the one 
of the dust, as presented in Sec. 3.1. The total mass of 
solid (defined as the mass in grains from 0.81 /xm to 8 mm) 
is constrained by the best dust model of Sec. 4.2 (M so iid 
= 4.9 x 1(T 7 M Q = 0.164 M e ) and we let the total 
atomic+molecular gas mass M gas vary from a gas-poor disk 
(M gas /M solid =10- 2 ) to a gas-rich disk (M gas /Af solid =10 3 ). 
The parameters of the model are summarized in Tab. 7. 
One of the main parameters controlling the gas energy bal- 
ance is the amount of PAHs (/pah), which remains uncon- 
strained in the disk of HD 181327 because their abundance 
is too small to produce detectable PAH emission features in 
the IR. We adopted two "extreme" values for /pah: a low 
PAH abundance (/pah=10 -5 ) and a high PAH abundance 

(/PAH=0.1) 3 . 

The non-LTE line fluxes for the [O I], [C II], CO J=3 
— > 2, and CO J=2 — > 1 transitions are given in Fig. 11 and 
Appendix B. Among all the lines, only the CO J=3 — > 2 
and J=2 — > 1 line fluxes behave monotonically, increasing 
in flux with increasing disk mass. The reasons for the appar- 
ent erratic behaviour of the fine-structure lines come from 
the non-linear nature of several physico-chemical phenom- 
ena. For example, the [C II] flux depends on the abundance 
of C + and on the excitation of the first energy level. The C + 
abundance first increases with total gas mass until CO self- 
shielding against photodissociation is efficient enough such 
that a significant fraction of C + is converted into neutral C 
and CO. The presence of CO also changes the heating and 
cooling balance, since the CO absorption of IR stellar pho- 
tons becomes the main source of heating. The change from 
a disk with no CO self-shielding to a CO self-shielded disk 
manifests itself as a jump in CO fluxes between a disk with 
gas-to-solid ratio of 1 to a disk with gas-to-solid ratio of 10. 
The abundance of PAHs has a weak influence on the gas 
temperature because the small UV flux from a F5/F6 star 
without excess accretion luminosity cannot ionize PAHs. 
On the other extreme, when the gas mass is very small (at 
gas-to-solid mass ratio below 1), the PAH and molecular 
absorption heating are nonexistent. Photoelectric heating 
by large grains dominates and explains why the line fluxes 
are similar for the high and low PAH abundance models. 

Line fluxes from models with /pah=10 -5 are consistent 
with the PACS upper limits for all gas-to-solid mass ratios 
(upper part of Table B.l). Only in the case /pah=0.1, the 
disk gas mass is constrained by [01] at 63 (im to be below 
~17 M e (M gas /Af solid < 10 2 , lower part of Tab. B.l). 

Eventually, even if the integration of the [O I] 63 and 
145 /im, and [C II] lines we performed was not sufficient 
to derive strong constraints on the gas-to-solid mass ratio 
in the disk, the combination of the GRaTcR and ProDiMo 



2 Woitke et al. (2009) provide a detailed description of the 
code, with subsequent improvements described in Kamp et al. 
(2010), Thi et al. (2011), and Woitke et al. (2011). 

3 The high PAH abundance is an upper-limit to the amount 
of PAHs in HerbigAe disks, where it is found to be of the order 
of 0.1-0.01 (e.g., Meeus et al. 2010, /pah=0.03 for HD 169142). 
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codes proved efficient to make predictions on the fluxes in 
the CO J=3 — s- 2 and J=2 — > 1 lines. These lines will be 
detectable for instance by ALMA, finally allowing to assess 
the gas content of the debris disk. 

Table 7. Setup of ProDiMo for the gas disk model 



Parameter Value 

Elemental abundances ISM-like 

UV field 1 Habing field (a) 

Cosmic-ray ionization rate £ 1.7 x 10~ 17 s _1 < - b ' 

Non-thermal line width dv 0.15 km.s^ 1 

/p AH (c) fir 5 - fir 1 

(a) Dalgarno (2006), (b) Standard interstellar UV field, (c) /pah 
scales with the gas mass and is defined with respect to an inter- 
stellar abundance of 3 x 10 -7 . 




0.01 0.10 1.00 10.00 100.00 1000.00 

Gas-to-solid mass ratio 



Fig. 11. Expected gas line fluxes as a function of gas-to-solid 
mass ratio for either a low-PAH abundance (thin lines) or a 
high- PAH abundance (thick lines). The horizontal lines denote 
the PACS 3(J upper-limits and do not include the 30% flux cali- 
bration uncertainty. The gas-to-solid mass ratios assume a solid 
mass of 0.164 M®. 



6. Additional constraints on the disk 

6.1. Scattered light: evidence for irregular aggregates 

The surface density profile used in this study comes from 
a new reduction of the scattered light images of Schneider 
et al. (2006) obtained with HST/NICMOS (see Sec. 3.1). 
This ensures by definition that the shape of any synthetic 
scattered light surface brightness profile in this study per- 
fectly agrees with the observed profile at 1.1 /im. However, 
the consistency is limited to the shape, and the agreement 
in terms of absolute flux is not guaranteed. This is because 
our fitting approach in Sec. 4 consists of adjusting the disk 
mass to get a best fit to the SED, which is equivalent to 
scaling up or down the surface density profile. A reason 
for adopting this strategy comes from the fact that the 
scattered light images do not directly provide the surface 
density S(r), but its product by the mean scattering cross 
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section er sca at the observing wavelength (1.1 /jm). The fit- 
ting process to the SED, on the other hand, yields an in- 
dependent, absolute scaling of the surface density profile, 
thus allowing one to obtain an initial estimate of <7 sca . We 
find (T sca = 3.7 /im 2 in the NICMOS 1.1 /im band. 

A second estimate of cr SC a is obtained by calcu- 
lating the theoretical scattering cross section for com- 
pact spherical grains given the grain properties (com- 
position, size distribution) that were inferred in Sec. 4: 
fTsca(A) = JaZaltin ™ 2 <9sca(a, A)dn(a) where Q sca is the 
scattering efficiency at wavelength A and for a grain radius 
a. For the best model (in the Bayesian sense, Tab. 6) we 
find <7 sca = 16.8 /jm 2 at 1.1 /im, a factor of about 4.5 higher 
than the first estimate. This implies a peak flux density of 
3.14mJy.arcsec _1 inconsistent with the 0.76 mJy.arcsec -1 
measured in the NICMOS profile. In other words, this 
means that our model predicts an albedo (w = 59% at 
1.1 /im) that is a factor of about 4.5 larger than the ob- 
served value. This suggests that smaller grains than those 
required to fit the SED dominate the scattering behavior 
of the disk. It is worth noting that such a conclusion is not 
specific to HD 181327, but is also observed in other debris 
disks. Recent examples include HD 207129 for which Krist 
et al. (2010) predict an albedo (based on dust properties 
also obtained from a fit to the SED) that is an order of 
magnitude larger than the value suggested by the scattered 
light image. A similar inconsistency is noticed for HD 92945 
by Golimowski et al. (2011) who point out that different ob- 
serving diagnostics suggest different minimum grain size. 

It is important to keep in mind that although our model 
suggests the grains are very porous and hence likely to have 
an irregular, fractal structure, we used hard spheres (Mic 
theory) to calculate both the thermal emission and scat- 
tered light optical properties. Unlike the thermal behav- 
ior, light scattering will likely strongly depend on the exact 
shape of the dust aggregates, and on the size of their el- 
emental constituents, two features that are not accounted 
for by the Mie theory (Voshchinnikov et al. 2007). Another 
piece of evidence that Mie theory might be a limiting fac- 
tor for a self-consistent model of HD 181327 comes from 
the observed scattering asymmetry parameter <7hg that is 
a measure of the direction of the light scattered by a par- 
ticle. For our best model (Sec. 4.2) the Mie theory yields a 
mean scattering asymmetry parameter |<7hg|= 0-95 at A = 
1.1 /im, much larger than the |<7hg|= 0.30 measured on the 
NICMOS image. Only a population of much smaller parti- 
cles could produce such a low scattering asymmetry param- 
eter if we stick to the Mic formalism, as already noted by 
Schneider et al. (2006). An analogous problem arises from 
the modeling of the debris disk of HD 92945 (Golimowski 
et al. 2011). Our micron-sized grains are highly anisotropic, 
essentially due to their spherical shapes. Incorporating vac- 
uum homogeneously inside the grains does not suffice to de- 
crease |<7hg| an d a more accurate description of the shape of 
fluffy aggregates is eventually needed. In fact, the scattering 
efficiencies and asymmetry parameters of fluffy aggregates 
are known to be correctly described by an EMT-Mie model 
(porous spheres) only in the Rayleigh regime, i.e. when the 
size of the inclusion (vacuum here) is small compared to 
the wavelength (see e.g. Voshchinnikov et al. 2007). This is 
not the case for the relatively large grain sizes and porosity 
when the wavelength becomes small. 



6.2. A narrower parent belt? 

Unlike scattering, we expect the continuum thermal emis- 
sion to be largely independent of the exact shape of the 
grains. Our best model of HD 181327's disk produces a 
sub-mm flux of 39 mJy at 870 /im, consistent within la 
with the LABOCA measurement (52 ± 6 mJy). The ATCA 
measurement at 3.2 mm was not model-fitted, because of 
large uncertainties on the optical constants at millimetre 
wavelengths. As a sanity check, we calculate the flux in the 
ATCA band using a power-law extrapolation of the SED. 
In the far-IR to mm regime, the model emission declines 
as A- 2 51±0164 . This results in a flux equal to 1.6t°; 4 mJy, 
which is within ~ 3 a compatible with the measured flux 
of 0.72 ± 0.25 mJy. Thus, the model remains valid for the 
longest wavelength observations (coldest emitting parti- 
cles), showing that the ring resolved in the near-IR with 
the HST is compatible with the spatial distribution of the 
largest grains, at least at first order. 

A limit of our approach is the assumption that the grain 
size distribution is the same everywhere in the disk, while 
it is known that only the smallest grains can populate the 
outer disk due to radiation pressure. The maximum size 
and power law slope of the dust size distribution should 
vary with the distance from the star (e.g. Fig. 4 in Augereau 
et al. 2001), whereas the model keeps both large grains in 
the outer disk and a constant slope for the size distribution. 
Such a radial variation of the size distribution is suggested 
by a comparison of a synthetic image at 70/im with the 
PACS 70/im map (Sec. 2.1). The synthetic image appears 
broader than the observed one with a FWHM of 7.68", 
versus 6.44" for the PACS image. Our interpretation is that 
the PACS band reveals larger grains that experience less 
the force due to the radiation pressure and thus are more 
bound to the parent belt. Assuming all grains are equally 
distributed results in putting too many large grains in the 
outer disk: a more complex model should account for the 
true dynamics of the particles. 

6.3. Emission from inside the belt? 

The region of the disk interior to 75 AU cannot be reliably 
assessed from the NICMOS data because of the corona- 
graph and artefacts from the PSF-substraction interior to 
1.5" (78 AU). Therefore we cannot exclude the existence 
of an additional dust component interior to the main belt, 
that would be invisible in the surface density profile we 
use. If micron-sized grains populate the inner disk, they 
would radiate mainly in the mid-IR, which might affect the 
results of our study. From the Gemini/T-ReCS images of 
Chen et al. (2008), we estimate that less than ~15% of the 
light emitted by the disk at 18.3 /im originates from a region 
interior to 1" (Fig. 2 in Chen et al. 2008). 

We perform a simple test of the robustness of our re- 
sults. Suppose that a belt of black body dust particles is 
orbiting the star at a distance r < 1". We calculate a black 
body spectrum than peaks at 18.3 /im (T = 150 K) and 
we adjust it such that its maximum flux represents 15% 
of the 53mJy attributed to the disk at that wavelength. 
Then we subtract the black body flux from the observed 
SED of Tab. 4 to simulate an observation of the system ex- 



4 Note that, at longer wavelengths, the emission seems to falls- 
off faster than a simple power-law. 
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eluding the hypothetical inner component. We reproduce 
the approach of Sec. 4.2 with this new SED. Interestingly 
the results of the Bayesian analysis are not significantly im- 
pacted. To account for the deficit of flux at mid-IR wave- 
lengths, only the minimum grain size needs to be slightly 
adjusted to eliminate the smallest grains: a min = 1.0 /im 
and k = —3.43, while we have a m i„ = 0.89 fj.ni and n = 
—3.42 for our best fit (Tab. 6). The composition remains 
unchanged for the best model (x? = 1.52) and no change is 
induced that would improve the consistency of the model 
with the scattered light flux (Sec. 6.1). Consequently we 
cannot reject the possible existence of an inner component 
to the disk, but none of our conclusions would suffer from 
its existence. 

7. Grain dynamics and ice origin 

7.1. Blowout size 

Grains in the HD 181327 debris disk for which radiation 
pressure overcomes gravity are put on unbound orbits and 
eventually ejected from the system on a short time scale 
(see for instance Krivov 2010, and ref. therein). The balance 
between radiation pressure (Fp) and stellar gravity (Fq) 
forces is evaluated with the distance independent, but size- 
dependent. /3 pr factor, given by 

a \Fpj_V (Qpr)£* 

Ppr ~ \F G \ ~ mAndGM* [ ) 

where (Q pr ) is the pressure efficiency averaged over the 
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Fig. 12. /3 pr — \Fp\/\Fg\ as a function of grain radius. The 
horizontal dashed line denotes the limit between bound (f3 < 
1/2) and unbound (/3 > 1/2) grains. 



stellar spectrum, m the grain mass, a the grain geometri- 
cal cross section, c the speed of light, G the gravitational 
constant, and and M* documented in Tab. 1. Assuming 
that the grains are spherical, Eq. 2 yields a critical size be- 
low which the grains originating from the parent belt are 
ejected from the system: abiow = a(/3 pr = 1/2) (assuming 
circular orbits for the parent bodies). 

In Fig. 12, we represent the /3 pr ratios as a function 
of grain radius for different compositions. Our best 4- 
component grain model has a blowout size abi ow = 4.83 /im, 
that is 5 times larger than the inferred minimum grain size 



in the system. Models implying no porosity nor carbon, on 
the other hand, have blowout sizes of the order of the in- 
ferred value of a m i n , more consistent with the theoretical 
argument that grains below the blowout size should be ex- 
pelled from the system. The response of silicate grains to 
the stellar irradiation and gravity fields is only slightly af- 
fected by the presence of ice, leading to similar shape and 
moderated values for /3 pr (a). Organic refractory grains on 
the other hand are highly sensitive to the radiation pres- 
sure and they can reach high /3 pr values. Fig. 12 shows that 
porosity has the effect of flattening the /3 pr (a) function, 
pushing the critical value (3 pl = 1/2 to larger grains. Hence 
the addition of porosity to a mixture of carbon and sili- 
cates produces a strong effect: /3 pr not only reaches a high 
value at its maximum due to carbon, but the f3 pi (a) func- 
tion decreases slowly such that the critical value f} pi = 1/2 
is reached only for large grain sizes. 

Decreasing the amount of carbon in the grains (it was 
imposed somehow arbitrarily, see Sec. 4.2) would mechani- 
cally decrease their blowout size. Pure silicate grains with 
the same porosity and ice content have abiow = 0.89 ^m, 
but they do not provide good fits to the SED. We com- 
puted additional grid of models with alternative vq/vsi 
values. The results confirm the conclusion of Sec. 4.1.1: all 
Si-r-C-r-H 2 0-|-porosity models lead to similarly good fits to 
the SED, as long as a sufficient fraction of carbon is con- 
sidered. For instance, if silicates represent only 10% of the 
C+Si mixture, then Vj co remains at 0.65 and the porosity 
is reduced to 55%, with Md us t = 0.043 M®, a m - m = 0.89 ^m 
and Xr = 1-66. All these values remain within the uncer- 
tainties of the original model. Models involving carbon-poor 
grains do not provide good fits to the SED, but as long as 
vc/vsi ^ 0.5, the results remain unchanged. In summary, 
the discrepancy between a m i n and abiow does not lie on the 
hypothesis we made on the composition of the grains. 

Our estimation of abiow is based on the assumption that 
the grains are spherical and homogeneous, but the actual 
response of inhomogencous aggregates to radiation pressure 
cannot be that simple. Saija ct al. (2003) for instance show 
that the fluffincss of the grains tends to decrease their ra- 
diation pressure cross-section, as if aggregates behaved like 
their elemental constituents. This would make grains below 
the theoretical blowout size less sensitive to radiation pres- 
sure, and allow them to remain bound to the star. Thus 
the apparent inconsistency may be interpreted as another 
signpost of the complex structure of the grains. We empha- 
size that fixing a m i n as equal to abi ow calculated with the 
(EMT-)Mie theory would not be a relevant approach when 
dealing with porous grains. 

Another possibility is that grains are not in dynamical 
equilibrium (i.e. an imbalance between production and ejec- 
tion), following for instance a major collisional event, such 
as an asteroid break-up. An excess of dust grains from all 
sizes would have been produced, including grains below the 
blowout size. These grains could have survived destructive 
collisions and radiation pressure blowing for long enough to 
be observed. Of course, observing such an event is unlikely, 
because dust typical lifetime is only of the order of 1,000 
years, which is small with comparison to the age of the star 
(see Sec. 7.2). But this "rare event hypothesis" would be 
consistent the scenario proposed by Chen et al. (2008) to 
justify the brightness asymmetry they claim to observe in 
thermal light. 
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7.2. Time scales and origin of the ice 

Collisions occurring between grains in a debris disk are a 
key process as they generally result in the production of 
smaller dust particles, either by fragmentation of the tar- 
get, or excavation. Any phenomenon longer than the char- 
acteristic collisional time scale will be considered negligible 
in the evolution of the dust population. Based on the prin- 
ciple that any particle will cross the disk midplane twice 
during one orbit, Thebault & Augereau (2007) proposed a 
simplified relation for the collisional time scale, valid for 
the smallest grains: t co \ = iorbit /f^ 60 , with i or bit the orbital 
period at the surface density peak, and rf^ the geometrical 
vertical optical depth. Here the orbital period at 89.5 AU is 
~ 850 years, leading to t co \ ~ 7600 years for the best model. 
This is more than three orders of magnitude smaller than 
the age of the system and two orders of magnitude smaller 
than the characteristic time-scale for Poynting-Robertson 
drag: t pT = 0.1 Myr for [3 W = 1/2 grains (Wyatt 2008). The 
evolution of the HD 181327 debris disk is clearly dominated 
by collisions requiring a mechanism to generate dust on a 
10 3 — 10 4 years characteristic time scale. 

Attempts to model time-dependent chemistry with 
ProDiMo reveals that the grain water ice mantle would 
be photo-evaporated in a few hundred years (for ~ 1 /jm 
grains), to a few thousand years (for ~ lOyum grains), a 
similar (or shorter) time scale to the dust mutual collision, 
requiring as well an ice reprocessing mechanism. This pro- 
vides a consistent view that the icy dust grains are main- 
tained through collisions in a reservoir of icy planetesimals. 
We highlight that, while detecting the mid- infrared spectral 
features in the emission of a debris disk is very challenging 
- unambiguous detection have only been made in brighter 
circumstellar disks (see, e.g. Malfait et al. 1999) - the sole 
analysis of the continuum SED alone can provide clues on 
the amount of ice in dust disks, thus on the nature of the 
parent bodies, often referred to as the leftovers of planet 
formation. 

In the Solar System Kuipcr Belt, the study of icy plan- 
etesimals, Kuiper Belt Objects (KBOs) in particular, is 
a prolific topic (for a recent review see e.g. Jewitt 2010). 
The presence of ice, together with the porous nature of 
the objects, is inferred from estimates of their densities. 
Unambiguous detection of water ice has been made possi- 
ble by the spectroscopic near-IR signature of crystalline ice, 
for instance by Delsanti et al. (2010) who completed their 
study of the plutino Orcus by radiative transfer consid- 
erations concluding that it contains a mixture dominated 
by amorphous and crystalline water ice. Confronted with 
space weathering, crystalline ice should be amorphidized in 
a short timescale. This implies that ice needs to be stored 
(and possibly produced) in the interior of the object before 
being brought to the surface, possibly through cryovolcanic 
events. The nature of water ice in KBOs remains unclear, 
but the vast majority of these objects are too small to ex- 
perience such geological processes (Dumas et al. 2011, and 
references therein). Thus KBOs covered by crystalline ice 
can only constitute particular examples, and the water ice 
covering the collisionally produced dust in the Solar System 
is more likely amorphous. The young age of HD 181327 al- 
lows the possibility that primordial crystalline ice has sur- 
vived in its debris disk. Our study tilts in favor of the 
amorphous nature of that ice, suggesting that either ice 
was predominantly incorporated in an amorphous form in 



HD 181327's Kuiper Belt objects, or that the amorphiza- 
tion timescale is shorter than ~10 Myr in this distant belt. 



8. Discussion 

The properties of the HD 181327 debris disk can be com- 
pared to those of the disks surrounding two of its co- 
moving "sister" stars, ft Pictoris (A3V, 1.7 M@) and 
AU Microscopii (M1V, 0.5 M@), which formed in the same 
molecular cloud from the same material and are coeval. 
Of course comparisons are limited by the specific models 
used to study each object. The relatively large porosity 
found here should be put in the context of the grain prop- 
erties found for AU Mic (Graham et al. 2007; Fitzgerald 
et al. 2007). Assuming a Si+C+H20 mixture, only highly 
porous grains (V = 80%), reproduces the scattered light 
profiles (exhibiting strong forward scattering) , the high po- 
larization and the thermal SED. In that model, the Mie 
theory remains valid to represent the scattering properties 
of porous grains, because the grains filling the outer disk 
are very small and stay in the Rayleigh regime. A total dust 
mass of O.OIMq is found, that essentially reflects the con- 
tribution from the inner mm-sized grains, concentrated in 
a ^40 AU belt. A similar porous grain model was used by 
Li & Greenberg (1998) to model /3 Pictoris yielding a dust 
mass ^0.4 Mg, contained in the ~ 100 AU-widc outer disk. 
Overall, "cold coagulation-like" dust models prove effective 
at reproducing observations of these young debris disks. A 
trend appears for increasing dust mass with spectral type 
as could be intuitively expected, but despite this difference 
of scale, the three stars are orbited by a similar material. 

Given the young age of these objects, we are investi- 
gating the conditions holding within disks just after the 
planetary embryos may have formed, the gas dispersed, 
and possible accretion onto outer giant planets happened. 
To explain the belt-like structure of the HD 181327 disk, 
the most straightforward explanation is to invoke the dy- 
namical influence of a planet inside the belt. The Wisdom 
(1980) analytical criteria gives the width of the chaotic zone 
a planet mean-motion resonances would create. Assuming 
circular orbits for the planetesimals, as suggested by the 
sharpness of the belt, this yields a range of possible val- 
ues for the planet semi-major axis depending on its mass. 
Assuming the inner edge of the disk is located at 78 AU 
from the star, the yet undiscovered planet would lie from 
73 AU for a Neptune-mass planet, down to 59 AU for a 10 
Jupiter-mass planet. 

The material at that stage appears as very hydrated 
providing a considerable reservoir of water in the outer 
system. The existence of such a reservoir is hypothesized 
in scenarios of the Solar System formation to explain the 
unexpectedly large water content on Earth for instance 
(Morbidelli et al. 2000). Indeed, terrestrial planets at this 
age (12 Myr) are likely too hot to hold on to water and 
later impacts of icy comet-like bodies from the outer re- 
gions arc a possible mechanism for delivering water and 
volatiles to planets in the habitable zone. Thus, not only 
HD 181327's debris disk reveals that the system has already 
grown planetary embryos - including maybe the bricks of 
rocky planets in the inner disk - but it likely holds the con- 
ditions that were prevailing before the Solar System Late 
Heavy Bombardment that depleted the Kuipcr Belt and 
determined terrestrial planets history. 
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9. Summary and conclusions 

Wc have presented new Hcrschcl/PACS and ATCA imaging 
of the HD 181327 debris disk, a 12Myr old star belonging 
to the j3 Pictoris moving group. The 70 and 100 /j,m images 
as well as the 3.2 millimetre observation reveal a marginally 
resolved extended emission. Fluxes in the three PACS pho- 
tometric channels were derived and complemented by con- 
tinuum fluxes obtained from PACS spectroscopic observa- 
tions. Together, these data provide precise information on 
the continuum thermal emission of the debris disk, but 
detailed models such as the one developed in this study 
are required to interpret the precise nature of the emit- 
ting material. Degeneracies between the grain properties 
and disk's spatial structure are broken in our model by 
making use of a reprocessed HST/NICMOS coronagraphic 
image revealing that the dust is bound to a narrow belt 
peaking at 89.5 AU from the central star. We focus our 
study on the dust properties (composition, size distribu- 
tion) using a radiative transfer code, and run a large grid 
of debris disk models for subsequent statistical (Baycsian) 
analysis. We conclude that the thermal light spectrum of 
the disk is well reproduced assuming a mixture of silicates, 
carbonaceous material and amorphous ice, as well as vac- 
uum to mimic a significant porosity. The observed SED 
can be very reliably reproduced by a population of grains 
larger than ~ 1 /xm, which follows a power-law distribution 
matching the ideal collisional equilibrium case and repre- 
sents ~0.05M® (1.5 x 10 _7 Af Q ) in grains with sizes up to 
1 mm. Extending the size distribution to the parent bodies 
of the observed dust grains leads to a mass of ~ 50M ffi for 
planetesimals up to 1 km. The composition we infer for the 
dust grains must be a good representation of a considerable 
reservoir of material hidden in unseen planetesimals. 

In this study, we attempted to use Hcrschcl/PACS spec- 
troscopy to detect for the first time far-infrared emission 
lines in the debris disk. The non-detections we obtain were 
used to investigate the gas content of the debris disk with 
our photochemical disk model ProDiMo based on the dust 
model found with GRaTer. The main uncertainty lies in the 
amount of PAHs in the disk. The current non-detections of 
[01] and [CII] lines alone do not provide unambiguous up- 
per limits to the gas content. Only in the case of a high PAH 
abundance, we can set an upper limit to the gas mass of 
~ 17 'M® (one Neptune mass). Coupling with other tracers, 
in particular with the CO lines accessible to ALMA, of- 
fers much better prospects to reach lower limits on the low 
gas content of debris disks. However, steady-state chem- 
istry is likely not ideally suited to debris disk modeling, 
and in future studies more appropriate prescriptions, time- 
dependent chemistry in particular, will be explored. 

We note that both the low albedo estimated from the 
scattered light images and the dust dynamics are not well 
modelled by spherical grains filled with vacuum, and that 
the grains are likely complex fluffy aggregates. Though their 
thermal behavior is correctly described by the EMT-Mic 
theory, the scattering of light at short wavelengths and 
their sensitivity to radiation pressure require a more accu- 
rate description of the exact shape of the aggregates. Our 
statistical analysis highlights that ignoring the porosity of 
the grains or setting the minimum grain size as equal to 
the blowout size would result in incorrect interpretation of 
the observations, in particular because the grain opacity is 
strongly dependent on the exact grain composition. 



We stress some necessary limitations of the models. 
First, we use a grain size distribution that is independent 
of the distance to the star, whereas it is known that only 
the smallest grains can populate the outer parts of the disk 
under the effect of radiation pressure. In this paper we dis- 
cussed clues that the observed scattered light profile may 
not be a perfect representation of that seen in thermal light, 
by Hcrschel or ATCA. The long wavelength emission em- 
anating from the largest grains might originate from an 
even narrower belt. Further studies will need to account for 
the radial differentiation of grain sizes with the objective 
to fit altogether resolved images from the near-IR to the 
millimetre domain. Secondly, the Dohnanyi-like power-law 
size distribution is known to be an insufficient description; 
if the real distribution implies wavy patterns as proposed 
by Thebault &: Augereau (2007), then adjusting a power- 
law likely puts too many /too few small grains in the disk. 
Dynamical models are required to understand better the 
general outcome of the collisional evolution of the disk. 
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Appendix A: Statistical analysis of the results 

The strength of the GRaTer code, in comparison with an 
optimization algorithm that would provide a unique best- 
fit solution and possibly miss local minima, is that all the 
models are stored enabling fast and easy post-processing of 
the results. 

To find the best solution in the grid of models, we use 
a classical least-square minimization approach with nf roc = 
3 to 4 free parameters depending on the grain composition 
and N = 48, which results in 43 to 45 degrees of free- 
dom (dof) (Tab. 6). The free parameters are k, o m i n and 
-Mdust, and depending on the situation vq, v- lcc or V. A/dust 
is handled separately because it is automatically adjusted 
to reproduce the surface density profile in the code. 

Let a be the vector of parameters. Here 
a = (k, a m j n , v lcc , V). For each set of parameters, 
we can define the quantity: 

a A (A - Mj) 2 

x w = 2^ 2 ( A - 1 ) 

t=i 0i 

where Di and M, ; are the observed and modeled values, 
respectively, and cr, is the measured uncertainty. For each 
X 2 value, we can define the likelihood of the data, i.e. the 
probability of the data given the parameters of a model. 
This assumes the errors are distributed normally and are 
uncorrelated. 

P(D|«) = ^«p(=^) (A.2) 
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where C is a normalization constant and D = {-D(i)}i=i...N- 
In order to evaluate the effect of each individual parame- 
ter a, on the quality of the fit, we implement a Bayesian 
inference method. We define the posterior probability of a 
as the probability distribution of the parameters a given 
the data D. It has to be noted that this formulation im- 
plicitly assumes that the model is suited to describe the 
observation. Baycs' theorem says that this quantity can be 
written: 



P(a\D) = 



P(a)P(D\a) 
P(D) 



(A.3) 



where P{D) is the prior probability of the data under all 
possible hypothesis (a normalization constant in practice), 
and P(a) is the prior probability distribution of the param- 
eters a. Here we use uniform priors for all the parameters, 
meaning that we do not favor any realization of the parame- 
ters prior to obtaining the data. That fundamental theorem 
provides a simple way to independently assess the param- 
eters of a model. Indeed, the marginal probability of each 
parameter a; can be computed regardless of the values of 
the other parameters a^: 

P{ ai \D) = J ... J J ... J P{ a] ^\D)da ...da^ 

(A.4) 

Then we obtain a probability law as a function of the re- 
alization of each individual parameter. Finally, we estimate 
the parameters by calculating the mathematical expecta- 
tion and variance. These are the values presented in Tab. 6 
where the uncertainties are 3<r. 



(A.5) 
(A.6) 



E{a { ) = J ai P( ai \D) 

ai min 

a 2 ( ai ) = E(a?) - E{ ai f 

Appendix B: Gas line fluxes 
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Table B.l. Gas line fluxes. The gas to solid mass ratios assume a solid mass of 0.164 Mj. The PACS 3-cr upper limits do not 
include the 30% flux calibration uncertainty. 
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